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The Green's Function for a Plane Contour. 

By Hobart Dickinson Frary. 



A Green's function is a function of a boundary and two points. Thei 
points, of which one is considered fixed and the other variable, both lie 
inside (or both outside) the boundary. In two dimensions the boundary 
is a contour. A large part of mathematical physics is devoted to the solution 
of boundary problems in two dimensions. All such problems can be solved 
immediately if the Green's function for the given boundary or contour is 
known. 

The character of the Green's function depends on three things, — (1) 
the differential equation of which it is a solution, (2) the contour or 
boundary to which it applies, and (3) the boundary conditions. The class 
of Green's functions to which we shall confine ourselves in this article is 
the one of most frequent occurrence in mathematical physics, namely, the 
Green's function for Laplace's differential equation, with what Hilbert calls 
boundary condition (I) ; that is, <7= on the boundary. Having found this 
particular solution G, we can solve the problem of Dirichet in the plane; 
that is, we can determine a solution u satisfying any given condition of the 
form u=U(6) on the boundary. This we can do by making use of the 
formula 

r dG 
u = 1/2tt ( U — ds, 

which in polar co-ordinates becomes 

/ 2V r dG dG 1 

There are two well recognized methods of obtaining the Green's function 
in two dimensions, both of which, however, are very limited in their applica- 
tions. These methods are the following: 

(1) The method of images. This method is easily applicable to the 
circle, semi-circle, infinite strip, half plane, equilateral triangle, 30 degree 
right-triangle, and the rectangle; but for other contours it would involve 
complicated Eiemann surfaces and hyper-elliptic functions. Apparently no 
one has attempted to apply the method to other polygons. 

(2) The method of Schwarz. Schwarz obtained a general formula 
w = /(«), where f(z) is in the form of a definite integral, by means of which 
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12 Fraby: The Green's Function for a Plane Contour. 

any given polygon can be mapped conformally on to the unit circle. Knowing 
this mapping function f(z), we can obtain the Green's function for the poly- 
gon immediately from the formula : 

= — Re log f(z) (Re = real part of) 

By this method we are able to get the Green's function for the equilateral 
triangle, the 45 degree right-triangle, an infinite strip, a regular five-pointed 
star, and a rectangle. Beyond this it is difficult to go. Other polygons lead 
to Abelian integrals and hyper-elliptic functions. Moreover, except for regu- 
lar polygons, it seems to be impossible to determine the unknown constants 
which appear under the integral sign in the function f(z). For regular 
polygons, however, we could evaluate Schwarz' integral, 

dz 

IV' 



/az 
„ (l — z n V 



(1 — *»)« 

in the form of an infinite series, and take the real part of the logarithm. In 
other words, the series obtained by expanding 

r z dz 
-? > 
0(1 — 2")" 

gives us the value of the Green's function for a regular polygon of n sides in a 
circle of radius 

dz 



f- 



o(l — 2")" 

where the point of discontinuity is at the center of the polygon. This work 
of Schwarz has been considerably discussed and amplified by other writers, 
particularly in a recent book by Study.* However, very little of fundamental 
importance has been added to Schwarz' original memoirs. 

The method of arriving at the Green's function which we shall employ 
in this article is new, although in principle similar to the methods used by 
Fourier and Neuman in attacking various physical problems. "We select an 
infinite set of linearly independent functions which are solutions of the given 
differential equation, and seek to expand the Green's function, or rather, the 
quantity — y, where y is the principal solution (log 1/r in this case), in a 
series of these functions, determining the coefficients of the series so that the 
boundary conditions are satisfied. The determination of these coefficients 
depends upon the solution of an infinite set of equations in an infinity of 
unknowns. In the problems solved by Fourier the solution of the infinite 
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set of equations was simple. In the problem before us the solution is much 
more involved, though entirely manageable in a large number of cases. 

This method of procedure is much more widely applicable than either 
of the two methods described above. These two methods give workable results 
for, at most, the regular polygons, two other polygons, the five-pointed star, 
and the rectangle. Our method is applicable, not only to all regular polygons, 
but to a large class of irregular polygons as well, and to many other contours. 
It gives the Green's function in the form, log 1/r plus an expansion in cir- 
cular harmonics, r n cos n6 and r n sin nO. The coefficients of this expansion 
are each given in the form of an infinite series of constants and can therefore 
be determined to any desired degree of accuracy by taking a sufficient number 
of terms. We are convinced that the method is actually workable in a practical 
sense in a large number of cases. 

We have confined our attention to a single class of Green's functions, but 
the same procedure might be employed to get the Green's function for any 
differential equation with any given boundary conditions.* It might also be 
extended to problems in three dimensions. 

Formulation of the Problem. 

By a plane contour we mean any closed plane curve which is regular in 
the sense of Osgood; that is, is composed of a finite number of analytic arcs 
or straight segments. An important case is that of the polygon. 

The Green's function for any such contour is given by the following 
formula in polar co-ordinates: 

Q(r,6)— logr + u(r,0), (1) 

where the function u is harmonic, that is, satisfies Laplace's equation, 

d 2 u 1 du 1 d 2 u 
dr^^r drr 2 W = ' 

and has no singularities in the region. 

Consider this function, u{r, 8). We know that it is everywhere harmonic 
in the area S enclosed by the contour G and has no singularities in 8. More- 
over, on the boundary, for which r = R(6) and u = U(6), we have G = and 

U{R,6)=\ogR{8). (2) 

Draw any circle, radius r (Fig. 1), about so that it lies wholly in S. By a 



* See Hadamard's Columbia lectures, pp. 47-52. 
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well-known theorem we can expand the harmonic function u uniquely in the 
series, 



u = AJ2 + 2 r n (A„ cos n6 -f- B n sin nO), 



(3) 



Where A n and B n are the Fourier constants for the area enclosed by the 
circle : ^ 6 ^ Sir, ^£ r 5= r . The functions r» cos n$ and r» sin nO form 
a complete set of orthogonal functions for the circular boundary. Let us write 
equation (3) for convenience in the form 



u = ^c m <j> m (r,6), 

m=o 

where tf> = 1, <j> ± = r cos 0, <f> 2 = r sin 6, 



(4) 




Fig. 1. 



Fig. 2. 



We can continue the function u analytically across the boundary C (Fig. 
1) by means of Schwarz' principle of symmetry. The extended function thus 
obtained takes values equal but opposite in sign at every pair of points which 
are symmetrical * with respect to an arc or segment of the. contour. At the 
symmetrical points or images 0' of the point 0, the function becomes infinite 
as log r, but at all other symmetrical points it is harmonic. Hence if we 
construct the images f of the point in the several arcs or segments of the 
contour, we know that the function u exists and is harmonic at least in a circle 
extending to the nearest image. Fig. 2 is drawn for a special case where the 
contour is a polygon. 

Now, the question arises: Is the function u represented by the same 



* For definition see Osgood, p. 671. 

f For the method see Lery's thesis on La Fonction de Green pour une Contour 
Algebrique, p. 49. 
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expansion in the whole circle ? To answer this question we cite the following 
theorem due to B6cher.* 

Theorem: If u is a function harmonic throughout the neighborhood of 
the point x y and if, when this function is continued analytically, the distance 
from x y to the nearest singular point of u, which lies in the same sheet of 
the Biemann surface generated by the analytic continuation in which x y lies, 
is K, then the development f in polar co-ordinates {the point x y being taken 
as pole), 

u = A /2 + S (A n r n cos n$ + B M r» sin n6), 

converges and represents u throughout the interior of the circle of radius 
r = K and does not converge throughout any continuum which does not lie 
in this circle. 

If we restrict the problem, then, to contours of such a nature that a point 
can be found whose images all lie outside of a circle about enclosing the 
contour (this we can make the unit circle), we have from the above theorem, 
for all contours of this type and all such points within each contour, the 
Green's function for the point and the given contour in the form 

00 

G = — log r + A /2 + S (A„r n cos n6 + B n r n sin n6) 
i 

00 

= — log r + ^c m <t> m (r, 6). (5) 



This expansion holds for all points inside the polygon and on the boundary. 
Moreover it is unique and we must be able to determine the constants A n and 
B n or c m from the boundary conditions. 

For the boundary we have r = B(6) and G = 0, that is, 

00 

A /2 + %[AJl»(6) cos n6 + B,Ji n (6) sinnfl] = log 5(0) 

i 

00 

or 3cm<M#) =logi8(0), (6) 

o 

where $ = 1 ; * 2 „-i = R n cos n6 ; * 2 „ == R n sin n$. 

Now we can apply a method of analytic continuation due to Borel, to 
obtain a new expression for the function u which will hold for a much more 
extended region, as follows: Equation (3) may be written 

u = A /2 + Re%{A n — iB n ) z n , 



* Trans. Am. Math. Soc, Vol. 10, p. 278. 

t BOeher gives here also an equivalent in terms of polynomials in ce and y. 
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where z is a complex variable. Since u is harmonic the expression 

00 

%{A n — iB n ) z n is analytic in a circle extending to the nearest image 0' 

i 

of 0. Consider for example the polygon* ABODE (Fig. 3) with the point 

of discontinuity at a point inside it. By means of Borel's integral we can 

00 

continue the function % A n z n analytically to the region enclosed by the polygon 




Fig. 3. 



A'B'G'D'E' whose sides are lines through the several images 0' and perpen- 
dicular respectively to the several lines 00'. The expression for the function 
thus extended is 



/(*)=/< 



00 00 






o n! 
Applying this to the problem in hand, we have, in the real plane, 

00 00 

u = A /2 + Ee C,-ty Un-iB n ) (zt)» 



n! 



t n 



= A Q /2 + fe"*Y (Ant* cos n6 + B n r" sin n6) — dt. 

1 

This function is harmonic inside the polygon A'B'G'D'E' whose sides are 
respectively parallel to the sides of ABODE. Hence this gives us a formula, 



00 00 



t» 



<7= — logr + A /2+ fe-'V (4„r» cos n6 + B n r n sin nO)—dt, (1) 
J i-mt n! 



Some or all of the sides of the polygon could be replaced by analytic arcs. 
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for the Green's function for any convex polygon whatever, with the point of 
discontinuity at any point whatever inside the polygon. Evidently this 
method can be applied, not only to polygons, hut to any convex contour what- 
ever and the point of discontinuity may be taken at any point whatever inside 
the contour. 

The whole problem is then reduced to the method of determining the 
infinite set of constants A n and B n , or c m . This problem may be solved in 
two ways: (1) by integral equations and (2) by an infinite set of algebraic 
equations. The first method gives only a formal solution; the second is a 
practical method as we shall show. 

Formal Solution by Integral Equations. 
The c's are Fourier constants. Hence, by a well-known theorem of 

00 

Fischer and Eiesz, the sum of their squares S c 2 m , converges. Moreover, if 

o 
we select any set of orthogonal functions whatever, linear, planar or solid, 

for example h m (x)> there exists one and only one function 2T(x) for which c m 

are the Fourier constants with respect to this orthogonal set. Whether the 

series thus obtained converges or actually represents the function at all is 

immaterial; for in any case the formulas for the coefficients hold, namely: 



H(x) Mx) d X> (8) 



where the definite integral is either single, double or triple, according as the 
orthogonal functions chosen are linear, planar or solid. Now substitute these 
constants back in equation (6) and we have, provided the series in the bracket 
converges, 

f #(x) [2 Mx) «*»(*) 1 dx=logR(0), (9) 

' A, L J 

an integral equation of the first kind with unsymmetric kernel K(x, 0) = 

00 

2Mx)*m(0). This integral equation has been studied by Picard, 
Schmidt, Bateman and others, and existence theorems have been worked out; 
but no general method of solution has yet appeared in the mathematical 
literature. The following special form of the equation, however, is of con- 
siderable interest: 

If we take for the set of orthogonal functions h m ( x ) the linear set cos n x , 
2 
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sinnx, the function E(x) is nothing else than the value of u on the unit 
circle, 

#(x)=m(i,x); 

and the kernel reduces to the well-known expansion for 

1 — B 2 

1 — 2J?eos( x — 9) + B 2 ' B<1 

Thus the integral equation (9) now becomes 

/2JT i B 2 (6) 

' u ^* h-m*~<x-!>)+*(') *-* m ' (10) 

where the left hand member is Poisson's integral. Hence we may state the 
solution of the problem before us as follows : 

Consider a closed contour of such a nature that a circle, with center at 
a point 0, can be drawn completely enclosing it, but not enclosing any of the 
images of with respect to the contour or arcs composing the contour. Make 
this the unit circle. Then the Green's function for this contour is obtained 
by using Poisson's integral backwards, as it were, to determine the value on 
the unit circle, namely H(x) = «(1, x)> of a harmonic function u(r, 6) which 
reduces to log B(6) for points of the given contour. Having found H we 
determine the Fourier constants or coefficients corresponding to it and substi- 
tute them in formula (5) to get the Green's function. 

On account of the complicated character of the kernel in equation (10) 
the solutioD of this equation seems to be impracticable. We shall therefore 
take up in the following pages an entirely different method for determining 
the constants for equation (5.) 

Practical Solution by an Infinite Set of Equations. 

Going back to equation (6), it is possible to reduce the problem of deter- 
mining the unknown constants c m to the solution of an infinite number of 
equations for this infinite set of unknowns. One way to do this would be to 
expand the various powers of B(6) as well as the function log B{6) in Fourier 
series and compare the coefficients of cos nO and sin nO on each side of the 
equation. This leads immediately to an infinite set of equations for the 
constants c m , but these equations are very awkward to handle. 

In order to get a more easily manageable set of equations we will con- 
struct a normalized, orthogonal set of functions from the set 

* = 1, $ 1 = .g(0) eos^, % = R(6) sinfl, % = B 2 (6) cos 26, , 
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functions which will be orthogonal for the variable 6 between and 2ir. 
These functions we can construct by the method of Goursat, provided the 
functions are linearly independent. We shall first investigate this matter of 
linear independence. 

An infinite set of functions *»(#) are said to be linearly independent if 
no finite number of them are linearly dependent; that is, if no equation of 
the form 

O o *r ,+ C x » ri + C 2 * r , + C 3 *r, + + C n *r= 0, (11) 

where the c's are constants not all zero, exists between a finite number of them. 
For the set of functions we are considering there certainly do exist such rela- 
tions for some contours r = R(6). For example, if B = sec 6 we have 

*o — *i-=0; 2% — $ 4 =0. 

However, for no closed finite contour can such a relation exist. For, suppose 
we consider a contour r = R{6) such that a linear relation of the form (11) 
exists. Then the left hand member of this equation is the value V(B, 6) of a 
harmonic function Y{r, 6) on a finite closed contour; hence by Gauss's mean 
value theorem, 

V(r,6)=0 

for all points inside the contour, which is only possible if all the coefficients 
vanish. Hence the set of functions 



$ = l, $ x = R eos 6,%=R sin 6, <J> 3 = R 2 cos 26, , 

and in fact any infinite set $»(0) obtained in the same way from a linearly 
independent set of two-dimensional harmonies <\> n {r, 6), are linearly inde- 
pendent for all closed contours, in particular for polygons. 

We can now apply Goursat's method and form a set, *<>> *i > 

of normal, orthogonal functions, as follows : 

% = a 20 $ + a 21 ^ x + a 22 $ 2 

* s = a 30 * + '<%*! + a S2 % + a 3S % (12) 



*m = Omo®o + tomi*i + «m 2 *2 + ««3*3 + + Omm$m 
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Solving for the *'s, we have 



*n 



1 = &„„*„ 



$! = B cos 6 <= 6 10 * + 6 u*i 

* 2 = B sin 6 = 6 20 * + b 21 % + 6 22 * 2 

$2«-l = Rn C °S ™0 = l 2 „_ 1} * + 6 2 «-l, !* x + & 2 »-l, 2 + (13) 

"T" 2 »-l, 2 »-l*2»-l 

* 2M = B n COS »10 = & 2 „, o*o + &2». 1*1 + 6 2 ». 2 + 

*T" 02». 2»-l*2»-l *T" ^2™, 2«*2« 



where 



6* m =(— l)"*" 1 - 



dm, m Mm+i, m+i ' 



' OJnn 



#»i+ 2 , m &m +2 , w+i Qm+2, m+2 " ' 



GSn-i, m ^w-i, w+i ^«-i, m+2 $«-i,m+ 3 ' 
&», m &», »»+i fln, >» +2 tt»,m+3 

thus : 







(14) 



Now expand log B(6) in terms of * , %, * 2 , 

logJB(fl) — 3U«*«, 

{=0 

where the coefficients are given by 

/27T 
logB(6)*i(6) d6. (15) 



Substituting in equation (6), we have 

c (&oo*o) + Ci (6 10 * + 611*1) 4-" c 2 (620*0 + 621*1 + 622*2) 

+ — 3Ui*i, (16) 

*=0 

where c , c 1} c 2 , are the unknown constants required. Comparing 

coefficients of the *'s on each side of the equation, we have 



"00 c o T "10 c i "T 620 C 2 -j- 6 30 c 3 -|- 

6n c x + 6 2 i c 2 + 6 31 c 3 -f- 

"22 ^2 "r "33 *3 "T 






(17) 



an infinite set of equations for the c's. Solving this set of equations formally 
for each c by eliminating the other c's one at a time, we have an expression 
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for each c in the form of an infinite series involving the 6's and the A's. If, 
then, we substitute for the &'s from equation (14), we arrive at the simple 
result 

c o == a oo-^o T a x<y°-x "T ^20-^2 T" a so-"-3 "T a 40-"-4 "T 

c x = a lx A t + a 21 A 2 + a 31 A 3 + a 41 A 4 + 

C 2 == *22-a-2 T a 32-^-3 ~T a 42^M "T 

(18) 

00 

Cm = ^j GimAi 
i=m 



Substituting this solution in equation (6), we have 

00 00 00 00 00 

S di Ai + 5 (5 «i, 2n . 1 AiB n cos n6 + S a«, 2 „AiB n sin «0) = 2 4i<Ev 

i=0 «=1 {=2»-l i=2» i=0 

It is sufficient to multiply equations (12) by A , A 1} A 2 , respectively, 

and add, to verify formally the correctness of this equation. Now since we 
know the c's exist and are unique, if we can prove that equations (18) give 
an actual solution of (17), this is the only solution. To prove this we shall 
make use of the following theorem due to Anna J. Pell : * 

Theorem: If the sequences {kit} and {j**} are such that the sequence 
{bih/kh} is of finite norm for every k, and the matrix (Mau/in) is limited, 
then for every sequence {An} such that {ju*A fc } is of finite norm, the system 
of equations 

00 

2 b ih d = A h (fc = 0, 1, 2, 3, ) 

has a solution C& such that {XfcC&} is of finite norm, and the solution is given by 

00 

c k = 2,a ih Ai (fc = 0, 1, 2, 3, ) 

Making A& = ju.i = 1, we can show that the three conditions required by 
the above theorem are satisfied, as follows: 

(1) The first condition to be satisfied is that the sequence {b ilc } shall be 
of finite norm for every h; that is, the coefficients in equations (13), taken by 
columns, must give sequences of finite norm. This we can prove as follows : 
Since the *'s are normal, orthogonal functions we have, for every n, 



* Ann. of Math., Vol. 28, p. 35. We have changed Mrs. Pell's notation to conform 
to our own. 
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S* cos n6 dO = ¥ 2n _ ±i „ + ^ ,+ + 6^, ^ (19) 

/27T 
A .B 2 »sinra0<Z0 = & 2 „ + 6 2 , + + 6 2 , 
2», ' 2», 1 ' ' 2)1, 2n . 

Hence for every fe, we have 

J27T /»2JT 

A i? 2 » cos «0 dO < I ..S 2 " <W = 2tr B 2n 
»/ • mam, 

J21T /»27T 

i? 2 " sin nO dO < I £ 2 » «W = 2ir 22 2re : 
o «/ o mas, ' 

so that we have, for the norm of the sequence {&«}, 

00 00 

SZ> 2 <4*S.B 2 " fc = 0, 1, 2, (20) 

But the maximum value of R(0) is always less than unity. Hence the right 
hand member of equation (20) is an absolutely convergent series. Thus we 
have proved that {&»&} is of finite norm for every Tc. 

(2) The second condition is that the matrix (au) of the coefficients in 
equations (18) shall be limited. This we can show to be true as follows: 
Consider first the matrix (&«) of equations (13). From (19) we have 

2 S & 2 < 4*r 2 B 2n < M, a constant: 

i=0 fc=0 ** 2„=fc •»«» 

that is, the &'s are of finite norm and hence * the matrix (&«) is limited. 
Now the matrix (a«) of equations (12) is the unique reciprocal of the matrix 
(bnc), and hence is itself limited.f Therefore its conjugate, the matrix (au) 
of equations (17) is limited. 

(3) The third condition is evidently satisfied; for the sequence {At} of 
the right hand members of equation (17) is of finite norm since the A's are 
the Fourier constants for the expansion of log.R(0) in terms of the *'s. 

Now by Mrs. Pell's theorem it follows that equations (18) give an actual 
solution of equation (17), and hence give us the required expressions for the 
constants c. 

Formula (5) for the Green's function may now be written 



G = — log r + A xf, + A^ + A 2 f 2 + A^ 3 + 
where 

'/'o = «00 

^i = ^io + a i\ r cos 



*Hellinger and Toeplitz: Math. Ann., Vol. 69, p. 307. 
t Hellinger and Toeplitz: Math. Ann., I. c, p. 311. 
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^2 = fl 2o + <h\>' cos 6 + a 22 r sin 6 

$3 = ft 30 + a 3i r C0S # + a 32 r S l n ^ + ^S** 2 C0S 26 



Now substituting for the coefficients A% from the equation (15), we have 

00 x»27T 

(? = — logr + 2^(r,0) j log 5(0) *i (6)^, (21) 

4=0 «/ o 

which, since the matrix (a«) is limited, can be written 

r27T 00 

logiJ(0 5 |>«(r, 0) 9i(t)] At. (22) 

Thus we have an expression for the Green's function in terms of a line integral 
around the "boundary. When r = R(0),ij/i(r, 0) = ** (6) , and the last mem- 
ber of the equation (21) becomes the expansion for log R in terms of the *'s, 
so that G vanishes identically. 

Application to a Square Contour. 

As a simple application of the above, let us find the Green's function for 
a square with discontinuity at the center. If we make the apothem of the 
square unity and measure angles from it, we have for the function R 

R(6) =sec0 — ,r/4= = „-/4 

J?(<9) =csc(9 „-/4 = =3ir/4 (23) 

R (0) = _ sec 6 3ir/4 = = 5ir/4 

R (6) = — esc 6 5^/4 = = 7,r/4. 

Beginning with the first function % = 1 we shall first orthogonalize all the 
other functions to it 'by determining C 2 n. t and <7 2 » to satisfy the following 
equations : 

J27T 
o (* 2 »-i — CVi^o) *o d9 = (24) 

J27T 
(* 2 „ — C 2 „* )* ^=0. 

From the first of these equations we have C 2n .i = f or n odd and w/2 odd ; 
that is, 

C 1 = a 3 = 5 = C 9 = C n = C' 13 = (7 17 =0. 

For n/2 even 

C 2 »-i = 2/ir. j see »0 cos n Q aft. 
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From the second of equations (22) we have (7 2 „ = for every n. Thus we 
obtain a new set of functions $', all of which are orthogonal to <5 , as follows : 

*,' = B cos $,' = B sin 

&, = B 2 cos 20 $ 4 ' = B 2 sin 20 

$ 5 ' = B 3 cos 36 $'„ = B a sin 30 

$; = i? 4 cos 40 — .48826 * 8 ' = i? 4 sin 40 

$ 9 ' = B* cos 50 ^ = i2 5 sin 50 

*,' = £* cos 80 — 1.11259 *' = B 8 sin 80 



*£ = i?" cos 120 — 3.01549 *' M = i? 12 sin 120 

Now we take *' and orthogonalize *,, *,, $ 4 ', to it in the same way 

as before. This process is continued until a sufficient number of orthogonal 
functions has been obtained. Finally these orthogonal functions must be 

r r 2K ~i J 

normalized by dividing each function * by I I <J> 2 d6 . Denote by * , * x , 

# 2 , , as in the preceding section, the normal, orthogonal functions 

obtained in this way. 

Next we determine the A's by means of the formula 



/27T 
log B(6)* m (6)d6. 




Finally, substituting in equations (18), we have the required constants for 
the square, namely, 

00 

c = 2 a i0 Ai = — .23653 —.05516 +.00883 +.00474 + • • • = —.28 • • • 
Ci = c 2 = c 3 = c t = c 5 = c 6 = 

oo 

c 7 = S a i7 Ai = —.11296 —.02100 —.00834 + = —.14 • • • 

c s = c 9 = c 10 = c u = c 12 = o 13 <= c 14 = 

00 

c 15 = S a ils Ai = —.00794 —.00013 -) = —.01 • • • 

4=15 

C 16 = C 17 = C 18 = C 19 = C 20 == ^21 B== ^22 ^ 

00 

c 23 = S a i23 A { = +.00153 + = +.002 • • • 

i=23 

Putting these constants in equation (5), we have the following formula for 
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the Green's function for a square with the point of discontinuity at the center 
and diagonal of length 2: 

Q = — \ogr —.28 —.14 r 4 cos 40 —.01 r 8 cos 80 +.002 r 12 cos 126 + ■ • • • 
= Re (— log z —.28 —.14 z* —.01 z* +.002 z 12 + ) (25) 

approximately. 

Thus we see that the expansion for u, that is, G + log r is the real part 
of an expansion in fourth powers of z. This may be verified by examining 
the expressions in elliptic functions, obtained by the method of images, for 
the square with the point of discontinuity and the origin of co-ordinates at the 
center. The character of the expansion for the Green's function for the square 
suggests immediately that for a regular polygon of p sides (p any positive in- 
teger) with discontinuity at the center, u will be the real part of an expansion 
in pth powers. It is evident at once that this is true when p becomes infinite, 
for in that case we have a circle for which u is constant. This may be shown 
to be true in general by making use of Schwarz' mapping function for a 
regular polygon of p sides. 



